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Abstract 

Abelian orbifolds of C 3 are known to be encoded by hexagonal brane tilings. To date it is 
not known how to count all such orbifolds. We fill this gap by employing number theoretic 
techniques from crystallography and by making use of Polya's Enumeration Theorem. The 
results turn out to be beautifully encoded in terms of partition functions and Dirichlet series. 
The same methods apply to counting orbifolds of any toric non-compact Calabi-Yau singu- 
larity. As additional examples, we count the orbifolds of the conifold, of the L aha theories, 
and of C 4 . 
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1 Introduction 

Brane tilings fH EJ have met with a lot of interest in the past few years. Each brane tiling 
gives rise to a quiver gauge theory, which can describe either the theory living on D3 branes 
probing a toric Calabi-Yau-three singularity or the theory living on M2 branes probing a 
toric Calabi-Yau-four singularity |3|. Faces, edges and nodes of the brane tiling - a periodic 
bipartite tiling of the plane - correspond respectively to gauge groups, chiral bifundamental 
fields and interaction terms in the superpotential. A periodic quiver can be constructed from 
the brane tiling by substituting nodes by faces and edges by arrows. The faces of the periodic 
quiver thus represent terms in the superpotential with a (negative) positive sign for (anti-) 
clockwise orientation and are in fact extensions of the usual notion of a quiver which does 
not have a "built-in" superpotential. 

Recently, the question of enumerating all possible brane tilings has been raised [4j. A 
classification of all brane tilings with up to and including Nj = 6 terms in the superpo- 
tential was given in this paper, and the results for all brane tilings with Nj = 8 terms were 
computed but not published. These results were derived using a computer code that reaches 
its limits for Nj = 10 terms. Thus a need for a better algorithm or a different approach arises. 
One possible approach is to count the number of tilings for a fixed number of terms Nj and 
collect the answer into a generating function. This turns out to be a difficult task and to 
date an answer is still unknown. We can therefore simplify the problem further and at- 
tempt to count the number of "sub-tilings". For example, we can ask how many tilings with 
n hexagons in the fundamental domain there are. Or we can ask how many inequivalent 
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tilings with 2n squares one can construct. Such questions turn out to be relatively easy to 
solve and are the subject of the present work. 

It is known that all orbifolds of a given geometry correspond to a repetition of the fun- 
damental domain [1J. For example, the number of inequivalent hexagonal tilings with n 
tiles equals the number of inequivalent orbifolds of C 3 with an Abelian group of order n. 
Similarly, the number of tilings with 2n squares is the number of inequivalent orbifolds of 
the conifold by an Abelian group of order n, etc. Note that in the case of compact Calabi- 
Yau manifolds, the problem is different because of the gluing conditions for the patches that 
result in a finite number of admissible orbifolds. 

It is useful to map the problem of counting brane tilings to the problem of counting sub- 
lattices. Take for example the problem of counting the Abelian orbifolds of C 3 . Instead of 
counting the lattices obtained by the repetition of n tiles, we can take the standard bipartite 
hexagonal lattice and count its sublattices of index n. Similarly, the problem of counting 
Abelian orbifolds of the conifold is equivalent to counting a certain type of square sublat- 
tices. We are thus led to the subject of enumeration of sublattices which has been studied by 
the crystallography community in great detail (see e.g. HHJIH I3IBI and references therein). 
Fortunately, we are able to take results from this field and apply them to the questions of 
interest in this paper. 

In the following we will outline some methods for counting the sublattices of a given 
lattice. While it is possible to enumerate by hand the first few orbifolds, this quickly becomes 
cumbersome. A general understanding of how the number of sublattices /(n) with index (i.e. 
size of the fundamental cell) n behaves and grows is therefore desirable. It turns out that 
this number decomposes into the symmetries of a given lattice. Our main results are closed 
formulae for the number of Abelian orbifolds of C 3 , of the conifold, of L aba theories and of 
C 4 , which are furthermore generalizable to any toric non-compact Calabi-Yau. In all the 
cases corresponding to Calabi-Yau-three geometries we find that for large n the dominant 
contribution is 

/(n) ~lcf' forn » 1 - ( L1 ) 

where G is the symmetry group of the fundamental cell of the brane tiling, and a is the 
sum-of-divisors function. 

The plan of this note is as follows. In Section [2j we outline the problem we are studying. 
In the following sections, we introduce the knowledge necessary to count all the sublattices 
of a given lattice and show how to apply it to the case of the hexagonal lattice (Abelian 
orbifolds of C 3 ). In Section[3j we discuss the cycle index which captures the symmetries of a 
given lattice. The Hermite normal form is introduced in Section [4] In Section |5j the concepts 
of the Dirichlet convolution and the Dirichlet series are introduced. In Section [6] we discuss 
the examples of the square lattice (Abelian orbifolds of the conifold), the L aha theories, and 
the tetrahedral lattice (Abelian orbifolds of C 4 ) following the same steps as detailed in the 
previous sections for the case of the hexagonal lattice. In Appendix|Aj we briefly discuss the 
general lattice in d dimensions. 
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n 1 2 3 3 



brane tiling 




geometry C 3 C 2 /Z 2 x C C 2 /Z 3 x C C 3 /Z 3 

Table 1: Brane tilings for the first orbifolds of C 3 . 

2 What we are counting and how 

To get a better understanding of our counting problem, let us give a brief description of some 
details. A more exhaustive explanation will appear in a forthcoming publication [9 J which 
will describe several computer codes that were used in order to obtain some of the results 
used in this note. 

We start by looking at orbifolds of C 3 . For simplicity let us focus on Z„ orbifolds. More 
general Abelian orbifolds such as Z„ x Z m can be treated in a similar fashion. Let us de- 
note the coordinates of C 3 by {21,22,23 }, and the orbifold action by (fli,fl2/<*3) such that 
{ Z\, 22, 23 } ~ { 0J ax z\, cv a2 Z2, to a3 Z3 }, with co 11 = 1 and a\ + a 2 + a 3 = mod n. In this nota- 
tion, the problem is to find all triples (fli,«2/ fl 3) that give inequivalent orbifolds of C 3 . The 
first few cases are as follows. For n = 1, the orbifold group is trivial and there is only one 
case, C 3 . For n = 2, there is again one case, which is commonly denoted in the literature 
as C 2 /Z2 x C. For n = 3, there are two cases, C 2 /Z$ x C with orbifold action (1,2,0) and 
C 3 /Z 3 with orbifold action (1, 1, 1). For n = 4, there are three cases, C 2 /Z 4 x C, C 3 /Z 4 with 
orbifold action (1, 1,2) and C 3 /Z2 x Z2, etc. The brane tilings for the first examples can be 
found in Table § a count of the first 16 cases can be found in the last row, / A , of Table g 

An alternative way of formulating the problem is by looking at the toric diagrams of these 
orbifolds. Since the toric diagram of C 3 is a triangle of unit area, the toric diagram of an 
orbifold of C 3 by an Abelian group of order n is again a triangle but with an area which is n 
times larger. The problem of counting all inequivalent orbifolds of C 3 is therefore equivalent 
to the problem of finding all triangles with vertices on integral points and area n. Of course, 
since these are toric diagrams, two triangles which are related by a GL(2, Z) transformation 
are equivalent. This provides another approach to counting orbifolds. A systematic method 
of finding such triangles is algorithmically different from that of finding inequivalent orb- 
ifold actions and therefore provides an alternative approach to the counting problem. As 
stated above, more details can be found in |9|. 

Yet a third approach is to think of the brane tiling as forming a bipartite hexagonal lattice 
and the problem of finding inequivalent toric diagrams is mapped to the problem of finding 
its sublattices. This is going be the language used in the next sections. The first step is to 
identify the symmetries of each lattice, which is the topic of the next section. 

There is however a subtle point that needs to be emphasized. The brane tiling encodes 
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the same information as the quiver diagram plus the superpotential, and in certain cases, 
one toric geometry admits several possible quiver gauge theories. In this case we speak of 
different toric phases. The lattices of the brane tilings pertaining to the different toric phases 
are different, but they preserve the same symmetries, and for the enumeration of toric ge- 
ometries, the symmetries are what matters. Even though our construction works starting 
from any lattice (brane tiling), if we want to count all the resulting quiver gauge theories, 
we need to keep track of the different toric phases that can appear in the process of orbifold- 
ing a given geometry. Consider for example the case of the conifold (see Section |&T) . The 
brane tiling is a bipartite square lattice. Dividing by Z2, one obtains Fo. This variety has two 
different toric phases, one corresponds to a sublattice of the square lattice, the other one is 
represented by a square-octagon lattice. While the quiver gauge theory of the former lattice 
is captured by counting the orbifolds of the conifold, the gauge theories stemming from the 
latter have to be considered separately. The final consequence is that even if we use the lat- 
tice of the brane tiling (which corresponds directly to the gauge theory), our counting covers 
the geometries but in general not all the gauge theories or toric phases that can arise from 
the orbifolds of a given geometry. 



3 Symmetries and the cycle index 

In the following, we need a way to capture the symmetries of a given lattice. Let us label the 
vertices of the fundamental cell by the numbers { 1, . . . ,m }. We now want to describe the 
group of permutations G of the set X = { 1, . . . , m } which result in the same fundamental 
cell. The cycle index encodes this information IfTOl . For the admissible symmetries, the fact 
that we are considering bipartite lattices plays an important role, since the preservation of 
the coloring of the vertices results in an additional constraint. 

In our case, it is most convenient to express the permutations in cycle notation. The cycles 
of g G G are the orbits of the elements £ G X under g. For each group element g we start with 
£1 G X and write down its orbit in parentheses, (ei g(ei) g 2 (£i) ■■■ g kl (^i))f where g k {e\) = 
Ei. We continue to do the same with the next element that has not yet appeared in an orbit 
until we have exhausted all the elements of X. Each g G G can thus be expressed in terms of 
«jt disjoint cycles of length k; cycles of length one correspond to elements that are fixed under 
g. The type of g is given by the partition of m [1^2 a2 . . . I* 1 ], where m = dl\ + 2*2 + ... + /«/. 
The partition is represented by the expression 

Zg(xi Xl )=x?x?...x?. (3.1) 

Definition. The cycle index of G is obtained by summing the over all elements g G G and 
dividing by the number of elements \G\: 

1 1 

Z G (xi,...,Xi) = — 1 Z g (xi,...,xi) = — T ^]c(ai,...,a/)^ 1 • ■ -zf', (3.2) 

where c(«i, ...,«/) is the number of permutations of type [l a i2 (l2 . . . l K >], and the sum runs over all 
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g Oil «2 «3 C(«j) £ 



/\ 1 3 - - 1 xf 

2 

3 

/\ (12) 11-3 xix 2 

1 

2 

/\ (123) - - 1 2 x 3 



3 1 



Table 2: Cycle index for the symmetric group S3. Z$ 3 = \ (x\ + 3xix 2 + 2x3). 
partitions ofm. 



Some cases of interest are: 



1. The cyc/fc group C„ is the group of symmetries associated to a circular object where 
reflections are excluded: 

z (c m ) = -LL<p( d K /d > (3-3) 



m A\ 

d\m 



where q>{d) is the totient function. 



2. The dihedral group D m is the group of symmetries associated to a circular object where 
reflections are allowed: 



Z(D m ) = l -Z{C m ) + < 



'l (m-l)/2 ., . , , 

jXix 2 , if m is odd, 

l/? (m-2)/2 . m/Z\ -r 

| ^qx 2 +x 2 J/ if m is even. 



(3.4) 



3. The symmetric group S m is the group of all permutations of n symbols: 

1 _ 

ui+2u 2 +-+kj k =m 1 «=1 71 k=l 

A convenient recursion formula is given by 

-1 m 

Z(S m ) = £ x k Z(S m . k ), Z(S ) = 1. (3.6) 

The link between the cycle index and the number of sublattices of the lattice L is provided 
by Burnside's lemma. 
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Lemma. Let G be a group of permutations of the set X. The number N(G) of orbits of G is given 
by the average over G of the sizes of the fixed sets: 

N (G) = t^t L \ F g\ ■> F g = {xeX\g(x)=x}. (3.7) 



The number f L {n) of sublattices of index n can be understood as the number of orbits 
of the symmetry group G when acting on the set X n of sublattices of index n. According to 
the lemma, this can be written as the average of the number of elements in X n that are left 
invariant by the action of g G G: 

/ L ( G ) = i^T E /» = l(^X„ | g(x) =x}\. (3.8) 

M g£G 

Using the cycle decomposition introduced above we can rewrite this expression as a sum 
over the types of the elements g, indexed by partitions a: 

f L (n) = T^E c Wti(n)- (3-9) 

By comparison with Equation 0.2\ we see that we obtain a subsequence for each monomial 
in the cycle index Zq(x). 

Every group contains the identity which is represented by the partition [l m ]. The cor- 
responding number of invariant sublattices f^m(n) only depends on the dimension of the 
lattice d = dim[L] and is given by the formula in Equation ( |A.10[ ): 

/£(»)= E hk 2 z ---k d d z{. (3.10) 

Example (A). Consider the bipartite hexagonal lattice corresponding to the geometry o/C 3 . Because 
of the bipartiteness , the symmetry group is not as expected for a hexagon, but S3, which we will 
also denote by a triangle. Table^shows the cycle decomposition for the symmetric group S3: 

Z S3 = ~(x 3 1 +3x 1 x 2 + 2x 3 ) . (3.11) 

Table^gives the number of sublattices of index n for the bipartite hexagonal lattice for each of the 
monomials appearing in the expression above. 

/ A = ^(/|+3/^ 2 +2/^) (3.12) 
The numbers / A (n), n < 500 are represented in Figure^ where the prime numbers are emphasized. 
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n 


1 


2 


3 


4 


5 


6 


7 


8 


9 


10 


11 


12 


13 


14 


15 


16 


J T 3 


1 


3 


4 


7 


6 


12 


8 


15 


13 


18 


12 


28 


14 


24 


24 


31 


r 1 
f A 


1 
1 


1 




2 
1 


3 
1 


2 



2 



2 
2 


5 



3 
1 


2 



2 



6 
1 


2 
2 


2 



4 



7 
1 




1 


1 


2 


3 


2 


3 


3 


5 


4 


4 


3 


8 


4 


5 


6 


9 



Table 3: Number of sublattices of index n for the hexagonal lattice, classified by the cy- 
cles of the symmetric group S3. According to the cycle index decomposition, f A = 

l(f x: • 3/ v , v , • 2/,). 




Figure 1: Scatter plot of the sequence / for a hexagonal lattice. Prime numbers are empha- 
sized in red. The two lines correspond to n/6 and e 7 n log log n/6. 
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4 Hermite normal form 



Consider a lattice generated by the d vectors (yi, . . .,y d ). Any sublattice L' of L d is gener- 
ated by d vectors {x\, . . ., x d ) which can be written as 



*2 



a n y 1 

a 21 y 1 + a 22 y 2 



(4.1) 



x d = a dx yi + a d2 y 2 H h a dd y d , 

where the integer coefficients satisfy the conditions 

< cijj < an V; < i . 



(4.2) 



The above construction is the so-called Hermite normal form. The index n of the lattice is given 
by the product 



n = Y\ a ii 

i=l 



(4.3) 



Expressing the sublattices via Equation (4.1 1, the problem of counting the sublattices of a 
generic lattice turns into the problem of counting the number of matrices fly that satisfy the 
conditions Equation (|4.2| and Equation (|4.3|. 



Consider a two-dimensional lattice L 2 . The condition a\\a 22 = n can be satisfied by the 
choice a 22 = m and a\\ = n/ m, where m is a divisor of n. If we want to count the number 
of sublattices invariant under the symmetry x a , we need to enumerate the possible values 
of a 2 \. The constraint a 2 \ < a 22 introduces a dependence of the number of possible values 
of «2i/ #{ a 2i } = <§y? (#22)/ on a 22 . The total number of sublattices f^(n) is thus given by 
summing ( m ) over all the divisors of n: 



(4.4) 



m\n 



Repeating the same construction for a three-dimensional lattice L3, we find that 

m\\nmi\m\ 



(4.5) 



A similar decomposition into d — 1 nested sums appears for any lattice of dimension d. 
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5 Generating functions 
5.1 The Dirichlet convolution 

Looking at Figure [TJ one realizes that prime numbers play a special role. This is one of the 
clues that point to the fact that the sequences corresponding to the monomials of the cycle 
index have the property of being multiplicative |TTT| . A possible explanation for this may be 
linked to the observation that the orbifold group Z pg is isomorphic to the group Z p x Z ? for 
p, q primes. 

Definition. A sequence f is multiplicative if 

f(nm)—f(n)f(m), when (n,m) = 1 , (5.1) 
where in, m) denotes the greatest common divisor between n and m. 

It follows in particular that / is completely determined by its values for primes and their 
powers, since for any n we can use the factorization n = p^ 1 p^ . . . p a / , and 

f(n)=f(p?)f(p?)...f(p a /). (5.2) 

Multiplicative sequences form a group under the Dirichlet convolution. For our examples, 
it is convenient to use this property and decompose each of the sequences into products of 
other sequences that are easier to deal with. First, we need the following 

Definition. The Dirichlet convolution of two sequences f and g is the sequence h defined by 

f(n) = (g*h)(n) = £g(m)h(^), (5.3) 

m\n 

where the notation m\n means that the sum runs over all the divisors m ofn. 

One can prove that this convolution is commutative, f * g = g * f, and associative, / * 
(g*h) = (f * g) * h, and that the sequence Id defined by 

Id(n) = {1,0,0,...} (5.4) 

is the identity, / * Id = /. To each sequence / one can associate its inverse f~ x satisfying 

/*/-!=/-!*/ = Id . (5.5) 

The inverse can be evaluated recursively via 

/-'(") = "T^y E W- (5.6) 

J v ' d\n 
d<n 

We have observed above that in two dimensions, the number of invariant sub lattices can 
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be put in the form of Equation ( |4.4| , 

/£(») = Es£K"0- ( 5 - 7 ) 

m\n 

This is equivalent to the statement that the sequence f^l is the convolution of gj; with the 
unit u defined by 

u(n) = {1,1/1/...}- (5-8) 



Its inverse is the Mobius function defined by 
u(n) = < 



1 if n is square-free with an even number of distinct prime factors, 
— 1 if n is square-free with an odd number of distinct prime factors, (5.9) 
otherwise. 



It follows that if / = g * u, then g = }i*f. 

Example (A). In the case of the bipartite hexagonal lattice we find: 

1. The sequence f 3 = { 1, 3, 4, 7, 6, 12, 8, 15, ... } corresponds to the identity permutation x\ and 
x i . . 

it is given by Equation d3J0> with d = 2. It can also be written as the convolution 



/i=u*N, (5.10) 

where 

N(n) = {1,2,3,...} . (5.11) 

2. The sequence f^ X2 = { 1, 1, 2, 3, 2, 2, 2, 5, ... } can fee written as the convolution of a periodic 
sequence of period 4 and £/ze unit: 

= { 1/ 0, 1, 2, 1, 0, 1, 2, 1, . . . } * u . (5.12) 

gi\x 2 zs zn the convolution of a finite sequence and u: 

/£* 2 = {1,0, 1,2, 1,0, 1,2,1,... }*u = {1,-1,0,2} *u*u . (5.13) 

3. TTze Zflsf sequence f X3 = { 1, 0, 1, 1, 0, 0, 2, 0, ... } also has the form of the convolution of the 
unity with a periodic sequence of period 3: 

f£ = {1,-1, 0,1, -1,0, 1,-1,0,... }*u . (5.14) 

The periodic sequence is the (non-principal) Dirichlet character of modulus three (see iUTI/ ): 

gx 3 =X3,2{n) = {1,-1,0,1,-1,0,...} . (5.15) 
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Putting all together we find that the sequence f A can be written as 

1 1 

fA = 6 ( f 4 + 3f ^ X2 + 2/ ^) = 6 (N +3 { 1/0/ ~ 1,2 } * U +2?C3 ' 2) * U ' (5 ' 16) 
5.2 Dirichlet series and power series 

The information contained in a sequence / can be usefully encoded into a generating function 
(more commonly used as a partition function in the physics literature). In the following, we 
will use two types of generating functions: 

1. the formal power series (partition function) 

oo 

f(0 = E /(*)'"; ( 5 - 17 ) 

n=l 



2. the Dirichlet series 



f (*) = E ^ • ( 5 - 18 ) 



n=l 

The corresponding inverse transformations are given by 

1 I F (t) 



(5.20) 



Dirichlet series are appropriate in the case of multiplicative sequences. In particular, if / 
is multiplicative, the series can be expanded in terms of an infinite product over the primes, 
the Euler product: 

™-£^-n(i + ^ + *£ + ...)- 

This is consistent with the observation that a multiplicative sequence is determined by the 
values taken for powers of prime numbers. 

Let us now consider the sequences that appeared in the example above: 

1. For the identity Id: 

Id(s) = l, ld(0=f; (5.22) 

2. For the unit u we obtain Riemann's zeta function: 
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3. For N: 



00 n 1 _i_ f3 

NW=E i5: =«-!)- N( t )=( T 3 7y t_-l; (5.24) 

4. For the finite sequence { 1, — 1, 0, 2 }: 

{1,-1,0,2X5) =1-1 + 1, {1,-1,0,2} (0 = f - f 2 H- 2f 4 ; (5.25) 

5. For the Dirichlet character X3,2, the corresponding Dirichlet series is the so-called L- 

function 



n=l 



Both types of generating functions have a simple behavior under Dirichlet convolution. 
Let f,g and h be such that 

f = g*h. (5.27) 

The power series for h reads: 



f(o = E /(»)»" = E E s( m ) *0 *" = E E *(*) • ( 5 - 28 ) 

n=l n=\m\n k=\m=\ 

This can be expressed in two ways, using the generating function for g or for /z: 

00 00 
F(0 = E ^(»»)H(t M ) = E fc(*)G(**) • (5-29) 

m=l fc=l 

In particular, since all our sequences can be written as sums over divisors (or equivalently 
as Dirichlet convolutions with the unit), we will always write 

00 

F(t) = ^G(t k ). (5.30) 

fc=i 

It is also possible to write the power series for the inverse of the Dirichlet convolution as 
follows. Let 

00 

f(t)= E g{m)h(k)t mk , (5.31) 

k,m=l 

then 

00 00 

H(t) = E m t k =E m g(k) F(t k ) , (5.32) 

k=l m=\ 

where ]i is the Mobius function. 
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More directly, the Dirichlet series is decomposed as 

/(«) ^^g(m)h(^) 



F s 



E 

n=l 



n=l m\n 



fc=l m=l 



(5.33) 



Note that the Dirichlet series corresponding to a sequence can be also understood as the 
Laplace transform of a discrete measure. This explains why it exchanges convolution and 
pointwise products. 

Both generating functions can be seen as linear transformations, hence the decomposition 
in Equation ( |3.9| still holds. 

Example (A). Let us now apply these formulas to the terms in f A : 

1. The sequence f^ is decomposed as f^ = u * N, hence the corresponding power series is gen- 



erated by 



and the Dirichlet series reads 



2. The sequence f^ X2 = { 1, — 1, 0, 2 } * u * u gives 



F$(s)=C(s)f(s-l). 



(i-0 (i + * 2 ) 



1-2 



-s i ol-2s 



1, 



(5.34) 
(5.35) 



(5.36) 
(5.37) 



3. TTze sequence f X3 = X3,2 * u groes 



F^(s) = L(s,^ 2 )as). 



(1 + f) (1 - t 2 



1 - f3 



1, 



(5.38) 
(5.39) 



Collecting all the terms we find: 
1. For the power series: 



G A (t) = 1 -G^(t) + 1 -G^ 2 (t) + 1 -G^(t) 



(1-0 (1 + t 2 ) (1 - f3) 



1, (5.40) 



whence 



m—l 



(1 - t m ) (1 + f 2m ) (1 - £ 3m ) 



00 00 



y y (— )" 2 f m ("i+ 2 " 2+3 "3) 

m=l ni,n2/"3=0 
^(0,0,0) 



(5.41) 
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symmetry Dirichlet series G(s) power series G(t) 
x\ C(s-l) 



(1 - t) (1 - *2) 

x:x 2 (1-2- + 2-) ^ {1 _)^ +t2) -l 

Table 4: Generating functions for the sublattices of the bipartite hexagonal lattice organized 
by symmetry. The actual sequences correspond to F(s) = G(s)£(s) and F(£) = EfcLi G(i fc ). 

2. For the Dirichlet series: 

(5.42) 



F A (s) = ^ (C(s - 1) + 3 (l - 2- + 2 1 - 25 ) C(s) + 2L(s^ 3/2 



The generating functions for the hexagonal lattice are summarized in Table^ 
5.3 Asymptotic behavior 

The asymptotic behavior of a sequence can be derived by looking at the corresponding 
Dirichlet series. 

Theorem. Let F(s) be a Dirichlet series with non-negative coefficients that converges for K(s) > 
cc > 0, and suppose that F(s) is holomorphic in all points of the line 3?(s) = oc, except for s = oc. If 
for s — > a. + , the Dirichlet series behaves as 

F(s) ~ A(s) + — — ~~~ ~~r , (5.43) 

where m S N, and both A(s) and B(s) are holomorphic in s = a, then the partial sum of the 
coefficients is asymptotic to: 

t»n~^«loz«\N). (5.44) 

n=\ 



In order to apply this theorem, we can make use of the following facts: 

1. The Riemann zeta function £(s) is analytic everywhere, except for a simple pole at 
s = 1 with residue 1; 

2. The L-function L(s,x) is analytic everywhere, except for a simple pole at s = 1 if x is 
a principal character. 

Another useful fact is Robin's inequality for the a function (sum of the divisors): 

a(n) < e 7 nloglogn, n large, (5.45) 
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brane tiling 



geometry 




conifold 



Table 5: Brane tilings for the first orbifolds of the conifold. 



where 7 is Euler's constant. This is true for large n, where large means n > 5041, and if and 
only if Riemann's hypothesis is true [12J. 



Example (A). The rightmost pole of the Dirichlet series F A (s) in Equation ( 5.42| | is found for s = 2, 
has order 1 and its residue is £(2)/ 6. Using the above theorem we conclude that the partial sum of 
the terms in the sequence f A behaves asymptotically as 



2 - 

72 



(5.46) 



and the sequence itself grows asymptotically as f A (N) = 0(N). 
For large n, the leading term is £(s)£(s — l)/6, hence 



rA , . e^nloglogn , 
f (n) < ° & , n large. 



(5.47) 



6 Examples 

6.1 Orbifolds of the conifold 

Another simple geometry that lends itself to the counting of its orbifolds is the conifold. In 
terms of the dimer model description, it corresponds to the bipartite square lattice with a 
black and a white vertex in its unit cell. The brane tilings for the first orbifolds can be found 
in Table [5] Note that as mentioned earlier, we are here enumerating all the toric geometries 
stemming from orbifolds of the conifold and not all possible quiver gauge theories, since we 
are not taking into account the multiple toric phases. 

The brane tiling corresponding to the conifold is a bipartite square lattice with two faces 
in the unit cell. Its symmetry is given by Klein's Vierergruppe V = Z2 X Z2, with cycle 
index 

Z v = ~ (4+24x 2 + 4) ■ (6-1) 

Accordingly, we can decompose / n as in Table |6j The numbers f a (n), n < 500 are repre- 
sented in Figure [2] 
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12 



Table 6: Number of sublattices of index n for the square lattice, classified by the cycles of the 
Vierergruppe V = Z2 x Z2. 




Figure 2: Scatter plot of the sequence / for a square lattice. Prime numbers are given in red. 
The two lines correspond to n /4 and e 7 n log log n /4. 
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symmetry Dirichlet series G(s) generating function G(f) 

*1 1 J (l-f)(l-f 2 ) 1 

*?x 2 | (L(s^ 4/2 ) + (1 - 2- + 2 1 - 28 ) ^(s)) (i- t )(i + *) " 1 

Table 7: Generating functions for the Abelian orbifolds of the conifold organized by the 
symmetries of the brane tiling. 

The analysis of the sequences corresponding to the terms in Zy gives the following re- 
sults. 

1. For the monomial x\ we obtain the by now usual term 

/° = u * N . (6.2) 

This is the same as for the example of the hexagon since it depends only on the dimen- 
sion of the lattice. 

2. The monomial x\xi receives two contributions (this is different from the other cases 
we have encountered): 

f% X2 = \ u * Oaa + { 1, "1,0,2 } * u) . (6.3) 
The corresponding power series is given by 

^-- (i-o'i + P) - 1 - < 6 - 4 > 

3. The monomial x\ corresponds to 

/S = {l,l}*u*u. (6.5) 
and the corresponding power series is 

G g<'>~ (1 . 1 ( )~(i% - 1 - < 6 - 6 > 

The generating functions in this decomposition are collected in Table [7[ 

Collecting all the terms and summing them according to the coefficients of the cycle 
index, we find 



1. For the power series: 
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Table 8: Number of sub lattices of index n for the L aha theories (Z2 symmetry), 
whence 



F°(f) = E 



m=l 



(1 - f«) (1 - f4«) 



00 00 



y y f.m(tli+4tl 2 ) 

m=l J3i,n2=0 
7^(0,0) 



(6.8) 



2. For the Dirichlet series: 



F n ( S ) = J as) [a« - 1) + 2 (1 + 2- 2s ) as) + a^)] 



(6.9) 



The asymptotic behavior of the partial sum of the terms in the sequence f a is given by 



N 



(6.10) 



and the sequence itself grows asymptotically as / n (N) = O(N). By Robin's inequality the 
sequence f u is bounded for large n by 



, n , . e^nloglogn 



(6.11) 



6.2 Orbifolds of the L fl&fl theories 



In the case of the L aha theories with a ^ b, the brane tiling lattice has Z2 symmetry and the 
cycle index is given by 



(*l+*2) 



(6.12) 



The number of orbifolds /^(n) (the first terms are collected in Table [8]) can be decomposed 
into two contributions: 



1. the usual term corresponding to the identity 

/£ = u*N; 

2. the term corresponding to the reflection x?. 

/£ = {!, -l,0,2}*u*u . 



(6.13) 



(6.14) 



18 



The power series reads: 



f*(0 = E 

m—1 



The Dirichlet series is 



1 + t 



3m 



(l-t m )(l- t 4m ) 



00 oo 



m=l ni,n 2 =0 

mo) 



(6.15) 



F*(s) = -£(s) l-2" s + 2 



>l~2s 



(6.16) 



From here we can read the asymptotic behavior of the partial sum of the terms in the se- 
quence f^: 



N 

L 

n=l 



71" 



-N 



4 24 

By Robin's inequality, the sequence / Q is bounded for large n by 

e 7 n log log n 



(6.17) 



n large. 



(6.18) 



If a = b, the brane tiling acquires an extra symmetry and the relevant group is Z2 X Z2. 
Since this is the same symmetry as for the brane tiling of the conifold that has been described 
in the previous section, the formulae in Equation (6.8 1 and Equation ( |6.9| apply also to this 
case. 



6.3 Orbifolds of C 4 

For the last example we consider Abelian orbifolds of C 4 . The three-dimensional counter- 
part of the brane tiling was described in IIT3I . The bipartite lattice has S 4 symmetry, like a 
tetrahedral lattice (see Figure 3J. For this reason we will denote the counting function f^(n) 
by a tetrahedron, just as before we used a triangle for the S3 symmetry of the orbifolds of C 3 . 
The cycle index for S4 is 

Z S4 = — (x 4 + 6 x\x 2 + 3 x\ + 8 x lX3 + 6 x 4 ) . (6.19) 

The first terms of the subsequences/^(ft) are collected in Table]^} and the first 500 numbers 
of f (n) are represented in Figure Ul 



Following equation ( |4.5| l, we know that each of the subsequences can be written in the 
form of a double Dirichlet convolution where one of the factors is the unity u: 

1. The sequence / 4 is the one corresponding to the identity. This means that it is given 
by Equation ( 3.10| for a = 3. Equivalently, it can be written as the convolution 

f$ = u * N * N 2 . (6.20) 
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Figure 3: The tetrahedral lattice. Points are at the vertices and the center of a tetrahedron. 
Four such tetrahedra connected at their vertices fit into a cube. 




Figure 4: Scatter plot of the number of sublattices of index n for a tetrahedral lattice. Prime 
numbers are given in red. The line corresponds to n 2 /24. 
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/° 1 2 3 7 5 10 7 20 14 18 11 41 15 28 31 58 



Table 9: Number of sublattices of index n for the tetrahedral lattice, classified by the cycles 
of the symmetric group S4. 
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2. The sequence /'f can be written as 



f$ X2 = { 1/ -1/0,4} *u*u*N . (6.21) 



3. The sequence coincides with /; 



x\x 2 ' 



4. The sequence fj~ x can be written as a convolution with the non-principal Dirichlet 
character of modulus three that has already appeared before: 



f£ X3 = { 1, 0, -1, 0, 0, 0, 0, 0, 3 } * u * u * X3 ,2 • (6.22) 

5. The sequence is the convolution of the non-principal character of modulus four: 

f$ = {1,-1, 0,2} *u*u*^,2- (6-23) 

We already have the Dirichlet series corresponding to each of these terms and we can collect 
them by using Burnside's lemma. The final result is: 



F V s ) = Ta (f? (s ) +6F ^^ +3F I + 8F ^3(«)+6^(s)) = 



m % l) (Us - 2) + 9 (1 - 2- + 2 2 - 2s ) as)) 



+ ^(8(i-3- s + 3 1 - 2s )l(^2,s)+6(i-2- s +2 1 - 2s )lU 4/ 2,s)) . (6.24) 

The rightmost pole is at s = 3 and has order 1. This means that the partial sum of the terms 
in the sequence behaves like 

V f*(n) ~ ^ (2)g(3) N 3 - 0.0274N 3 . (6.25) 
3 x 24 

It is also possible to write the power series corresponding to each term as follows: 

CO 

L G $(n, (6-26) 

m=l 
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symmetry 



Dirichlet series G(s) 



generating function F(t) 
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X X X 3 
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as-ms-2) 



£ nm 2 t mn 

n,m=l 



n,m=l 



n,m—l 



00 

y ^n 2 +4m 2 _ j 
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00 

y ^n 2 +mn+7m 2 
n,m=—co 



Table 10: Generating functions for the Abelian orbifolds of C 4 , organized by the symmetries 
of the brane tiling. 



where 



Gj(0 = £ nm 2 t mn , 



n,m=l 



00 



n,m=l 



G x x x 3 (0 - 2 
1 ' 



E ' 



n 2 +4m 2 



n.m=— 00 



g2(0 = 



E ' 



n +mn+7m 



n,m=—co 



(6.27a) 
(6.27b) 

(6.27c) 

(6.27d) 
(6.27e) 
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A Generic lattice in d dimensions 



In this appendix we derive a general expression for the number f Ld (n) of sub lattices of index 
n for a generic lattice in d dimensions. Using the construction in Section |3J we want to count 
matrices of the form 

/ An \ 

#21 a 22 

(A.l) 

\#di a d2 ... a dd J 
where flnfl22 • • • «dd = n, and < fly < an, V; < i. 

• For d = 1, the only condition is flu = n, so there is exactly one sublattice for each 
choice of the index: 

f Ll (n) = l. (A.2) 

• For d = 2 we want to count the triangular matrices 



flu 

fl21 «22 



(A.3) 



such that flnfl22 = n and «2l = 0, 1, . . fl22 — 1. It is immediate to see that there are m 
such matrices for each choice of a divisor moin: 



f L i(n) = Y j m= £m/ Ll (m). 



(A.4) 



In d = 3 we count the matrices 



m\n mm 



fa n o\ 

fl2l fl22 
\«31 «32 «33 / 



(A.5) 



where fln«22 fl 33 = n / fl 2i = 0, 1, . . . , fl22 — 1 and both and 032 are in { 0, 1, . . . , 033 — 1 }. 
Let us write 



flu 



mi 



«22 



ni2 



«33 = m 2 , 



(A.6) 



then we have ^22^33 = tn-i/mj x = m\tn2 such matrices for any divisor wzi of n and 
any divisor m2 of m\. Hence 



f l3 {n) = m i m 2 = ^m/ L2 (m). 



(A.7) 



For general one has the condition flnfl22 • • • ^dd — n an d a total of ^22^33 • • ■ fl 
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matrices. Writing the coefficients as 

n mi m d _ 2 , A ON 

«11 = / «22 = / ••• a (d-l)(d-l) = / a dd = nid-l , (A.8) 

nti ni2 v m rf-i 

we obtain 

/ Lrf (") = E E • • • E m 1 m 2 ---m d _ 1 = J2 m f L ^ 1 ( m )- ( A -9) 

mi\nni2\mi nid-i\ m d-2 m \n 

The generating functions for these sequences are written as follows. From the observa- 
tion that there are 022033 • • • d^ 1 matrices of index n = flnfl22 • • • ^dd we can obtain directly 
the power series 

00 00 
F L « (t) = Y J f Ld {n)t n = £ k 2 k 2 3 ■ ■ ■ rf-tyik-kt _ (A.10) 
n=l k\,ki,...ka=l 

The Dirichlet series can be written by using the recursion relation 

f L i {n) = Y J mf 1 "- 1 (m) . (A.ll) 



Rewriting it as 



n n/m 

r e see that the sequence [f Ld (n)/ 



Ei/ L "W, (A.12) 



we see that the sequence [f Ld (n)/n] is the Dirichlet convolution of the sequences [1/n] and 



[£^)] = [i]*^x(„)]. (A.13) 
It follows that 

F L <<(s + l) = ^(s + l)F L "- 1 (s). (A.14) 
This recursion relation can be solved starting from 

F L i(s)=£(s), (A.15) 

and we find 

F L *(s) = - 1) • • • £(s - d + 1) , (A.16) 
or in the form of a Dirichlet convolution 

f Ld = u*N * N 2 * • • • * N d_1 , (A.17) 

where N k {n) = { l k ,2 k ,3 k , . . . }. 
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